Fermi Surface Evolution Across Multiple Charge Density Wave Transitions in ErTe3 
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The Fermi surface (FS) of ErTes is investigated using angle-resolved photoemission spectroscopy 
(ARPES). Low temperature measurements reveal two incommensurate charge density wave (CDW) 
gaps created by perpendicular FS nesting vectors. A large Ai = 175 meV gap arising from a CDW 
with c* — qcDwi ^ 0.70(0)c* is in good agreement with the expected value. A second, smaller 
A2 = 50 meV gap is due to a second CDW with a* — qcDW2 ^ 0.68(5)a*. The temperature 
dependence of the FS, the two gaps and possible interaction between the CDWs are examined. 



PACS numbers: 71.45.Lr, 79.60.-i, 71.18.+y 

Charge density wave systems have been studied for 
many decades [1 . CDWs and spin density waves, the 
spin analog, exist in all dimensions. However, enhanced 
electron interactions in systems with reduced dimen- 
sionality increases the susceptibility for CDW forma- 
tion, also known as a Peierls distortion in one dimen- 
sion [liaillSlElIIlE]. In ideal ID systems the FS 
topology determines the susceptibility for CDW forma- 
tion via electron-phonon coupling [1 . If the FS can be 
nested with one q- vector of a particular phonon mode, 
the ground state energy can be reduced by electron- 
phonon coupling resulting in gaps opening at the Fermi 
level. Real systems are never perfectly ID, resulting in 
imperfect FS nesting with a partially gapped FS and 
residual metallic pockets [6", "71 [8] . While a great wealth 
of information has been learned from CDW systems, mi- 
croscopic models explaining CDW mechanisms have been 
elusive and debated [H [TOl [TTl [12]. It has even been 
proposed that the idea of FS nesting is insufficient to 
explain the nesting vectors realized in the family of rare 
earth tritelluride (RTes) compounds [11 . The RTea fam- 
ily offers a unique opportunity to systematically study 
CDW formation over a wide range of tunable parameters 
via rare earth substitution [12l [131 HH [13 US]- ad- 
dition, the simple electronic structure makes RTes niore 
tractable for theoretical modeling to gain a deeper un- 
derstanding of the CDW phenomena and electron corre- 
lations in general [H HI [17] . 

The CDW properties for several members of the RTes 
family have been well characterized previously (T] [H [131 
[Ml [T5l [TB] . RTes has a layered, weakly orthorhom- 
bic crystal structure (space group Cmcm) consisting of 
two planar Te nets sandwiched between buckled RTe 
slabs. X-ray analysis reveals an incommensurate lat- 
tice modulation characterized by a single wave vector 
(g - 2/7c*) HI [E]. Previous ARPES studies show a 
partially gapped FS with the gap maximum occurring 
along the c-axis (i.e. V — Z or kx =0 since the 6-axis 
is the long axis) [71 El [161 [E] • The gap evolution as a 
function of R is consistent with a FS nesting driven si- 
nusoidal CDW with the bandwidth and density of states 




FIG. 1: ErTes FS determined from ARPES. Data taken at 
T = 10 K with hv — 23.5 eV and a 10 meV integration 
window around Ep- Solid square outlines 1** BZ of the true 
3D crystal structure. Dashed lines are expected FS from 2D 
Te net unmodulated by CDWs. Red dotted square outlines 
FS areas shown in Fig. [3] 



at the Fermi level {N{Ef)) tuned by a shrinking unit 
cell [I51[I1II71. 

The recent discovery of a second CDW transition in 
heavier members of RTes is a rare display of the riches 
in CDW formation and a unique opportunity to further 
test the concept of a FS driven CDW. The second CDW 
was ffist discovered by transport measurements [15 . It 
has further been conffimed by x-ray scattering displaying 
a second lattice modulation perpendicular to the ffist, 
parallel to the a-axis with q2 ~ l/3a* [15j. We have 
performed ARPES investigations of ErTes, specifically 
chosen since it has two well separated CDW transitions 
at 155 K and 267 K. Our finding further strengthens the 
notions of FS nesting driven CDW formation in this par- 
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ticular compound. Data suggest that as the rare earth 
ion is varied, the second CDW is formed only when the 
first CDW weakens with the decreasing lattice param- 
eter, making larger FS segments available for the new 
nesting condition to form ^15| il6| . While the result can 
be qualitatively explained by a simple tight binding (TB) 
model, the data contain rich subtleties suggesting a dy- 
namic interplay of the two CDWs. Our finding establish 
ErTea as an excellent model system to study the evolu- 
tion and entanglement of two many-body states existing 
within the same atomic plane. 

Single crystals used in this study were grown by slow 
cooling a binary melt and have been well character- 
ized elsewhere [15]. All ARPES data were taken at the 
Stanford Synchrotron Radiation Laboratory beamline 5- 
4 with an energy resolution of 10 meV and angular reso- 
lution of 0.3°. 

The FS determined from ARPES is shown in Fig. [l] 
A large gapped region in the vicinity of the T — Z axis is 
consistent with previous studies of RTea members with 
a single CDW. While x-ray data shows qcDW2 to be 
parallel to the a- axis, the FS along T — X {k^ = 0) is 
not completely gapped at T = 10 K as shown in Fig. [T] 
The intensity is weak along the inner FS, but the outer 
FS retains significant spectral weight. Further inspec- 
tion slightly off the T — X axis {k^ ~ 0.07c*) reveals a 
second gap does appear in both the inner and outer FS 
pieces. In addition, several weaker features are evident in 
the regions around the large and small gaps just below 
Ep. Fig. is of the same data as Fig. [l] but with a 
logarithmic intensity scale and slightly larger integration 
window to emphasize the weaker features near Ep. As 
shown later, the weaker features provide crucial evidence 
for the induced FS from both CDWs. Given the complex- 
ity, there is amazing degree of agreement between data 
and the expected FS pieces. 

While ARPES exposes distinct differences between the 
FS of lighter RTes members with a single CDW and 
the FS of ErTea [3 13 [16 , differences in the theoreti- 
cal FS are subtle [15]. Linear muffin tin orbital (LMTO) 
calculations suggest the electronic structure near Ep is 
still dominated by the Px and Pz orbitals of the two Te 
planes [151 Ull IE]- Hence, one would expect that the 
TB model of a single Te plane which so accurately de- 
scribes other rare earth members would also accurately 
describe ErTes. Details of the TB model are described 
elsewhere [8 , T6^. A square net of Te atoms is assumed 
with perpendicular chains created by the in-plane Px and 
Pz orbitals {py orbitals are assumed to be filled). The 
model parameters consist of electron hopping terms along 
a particular chain (ty) and perpendicular to the chain 
{t±). The Fermi velocity, the slope of the bands travers- 
ing Ef, is primarily determined by while a FS cur- 
vature is introduced proportional to t±/t\\. In this work 
the TB model is expanded to not only include effects of 
the \/2 X \/2 i?45° reduced 3D Brillouin zone on a 2D 
Te plane and ±(c* — qcDWi) as used previously, but in- 
cludes the second CDW with the addition of Px/Pz bands 




0.0 0.1 0.2 0.3 § -1.2 -0.8 -0.4 0.0 0.4 
kx/a* E - Ep (eV) 



FIG. 2: (a) ARPES FS (top) compared to interacting TB 
model FS (bottom) with either one (left) or two (right) 
CDWs. Both the data and model use a logarithmic intensity 
scale and 20 meV integration window about Ef to emphasize 
weaker features just below Ef- The TB model was convoluted 
with an angular resolution function. Vertical dotted line dis- 
tinguishes between TB model with one or two CDWs. (b) 
Energy of leading edge gap (LEG) in the dispersion relative 
to Ef measured along the inner FS for a constant kz- Note 
the CDW gap magnitude is measured relative to the center of 
the gap (see Fig.jsj). The markers are the data and the lines 
are from the TB model. Shaded areas correspond to regions 
where entire CDW gap is below Ef resulting in intensity at 
Ef (see text), (c) Density of states from the TB model with 
different number of CDWs present. 



translated by ±(a* — qcDW2)- 

To further understand the implications of two trans- 
verse CDWs on the FS, the TB model is further ex- 
tended to include interactions between various bands. An 
electron-phonon coupled Hamiltonian is used to model 
the band interactions in a perturbative fashion [TJ [16]. 
Interactions arising due to the 3D crystal structure and 
both CDWs are included in the Hamiltonian, however, 
only first order interactions due to band crossings arising 
from a single (^-vector translation are included. While 
higher order terms are expected to exist due to the 
incommensurate nature of the CDWs, the intensity is 
expected to be extremely weak and thus can be ne- 
glected [8l[l6l[20]. The lattice parameters {a = c = 4.27 A 
assuming a square Te net) are taken from x-ray results 
and all other model parameters are estimated from the 




data, both from the FS and band crossings occurring 
below Ef. Since the intensity appearing in the recon- 
structed FS originates from the original Px and bare 
bands, the color scale in the model FS is proportional to 
the sum of the square of the resulting Px and Pz eigen- 
vector amplitudes from the interacting Hamiltonian [16] . 
The resulting reconstructed FS from the interacting TB 
model is shown in Fig. [2^. The spectral intensity existing 
near the FS corners at both kx = and = are now 
explained by the electron-phonon coupled TB model of 
the CDW. Shadow bands, bare Px and Pz bands trans- 
lated by CDW g'-vectors, correspond to weak features 
observed in the data revealing the two transverse CDWs. 

Due to the agreement between the model FS and the 
observed one, the extended interacting TB model is used 
to characterize the CDW properties in a similar man- 
ner to previous studies of lighter RTea members [16]. 
The gap evolution along kx for the inner FS square is 
shown in Fig. [2]3, demonstrating the effects of multiple 
CDWs. The excellent agreement between the model and 
experimental data suggests the observed FS is consistent 
with a nesting driven CDW, however, quantitative dis- 
crepancies exist which must be discussed. The model 
intensity at the corners of the inner FS square do not 
match the experimental data. This discrepancy is ob- 
served both along the F — X and T — Z directions, but is 
most pronounced along T — Z. This discrepancy is most 
likely due to the simplicity of the model, neglecting ef- 
fects from higher order CDW terms, interactions beyond 
the nearest neighbor in-plane p-orbitals, a bi-layer split- 
ting arising from neighboring Te planes and orthorhom- 
bicity of the unit cell [16j. Despite these discrepancies, 
the agreement justifies using such a simple model for il- 
lustration purposes. The density of states is calculated 
for the interacting TB model and shown in Fig. The 
onset of a single CDW suppresses N{Ef) to ~ 77% of 
the unmodulated value while the second CDW further 
suppresses N{Ef) to ~ 74%. Although the gains due to 
the second CDW are modest, mean field transition tem- 
peratures depend exponentially on N{Ef)- In addition, 
estimations of the Te bandwidth for GdTea (c = 4.33A) 
is 4.70 eV [H] while for ErTes (c = 4.29A) it is 4.85 eV. 
While ErTes ARPES data is the first in the RTes family 
to show multiple gaps in the FS, the increasing band- 
width proportional to the shrinking lattice agrees with 
the general RTes trend. 

Studies of lighter RTes members with a single CDW 
suggest a shrinking gap accompanies the shrinking lattice 
parameter, allowing a larger portion of the FS to remain 
intact. In addition, the gap maximum always appears 
along the c-axis [16 . Perfectly nesting the corners of 
the FS along each axis (i.e. c* — qcDW 0.68c*) will 
pin the center of the CDW gap at Ef [E]. However, 
observed nesting vectors and hybridization of Px and Pz 
bands effectively pushes the center of the gap below Ef- 
This effect is most pronounced at the corners of the inner 
and outer FS allowing intensity to be observed in these 
regions. The gap magnitudes (Ai = 175 meV and A2 = 




(d)0.1 

> 

0^-0.1 



LU 



-0.2 



-0.3 



'iOK 


iSOK 


250K 


275K 


'ft 


i 


i 





(e)0.i 



Wave Vector 



■D 


N 





1 








J \ \—LJ L 



-0.4 



-0.2 -0.2 -0.1 

E-Ep (eV) 



50 100 150 200 250 300 

Temperature (K) 



FIG. 3: (a-c) ErTes Fermi surface taken at different tem- 
peratures. The data shows one section of the inner FS to 
emphasize the effects of the closing gaps. Markers correspond 
to locations of the measured gaps. Horizontal lines in (a) rep- 
resent extent of cuts shown in (d-e). (d-e) Spectra divided 
by the Fermi function at different temperatures for the large 
and small gaps respectively, (f-g) Temperature dependence 
of the fitted and denoised dispersion data for the large and 
small gaps respectively. Different colored lines correspond to 
different temperatures. Dashed vertical lines mark the center 
of the closing gaps, (h) Summary of the two energy gaps for 
different temperatures. Solid lines show expected mean- field 
behavior with energy gaps and transition temperatures scaled 
to experimental values. Note that the magnitude of each gap 
is measured from the center of each gap marked by horizontal 
dashed lines. 



50 meV) are smaller in ErTes than the lighter rare earth 
family members (for CeTea A ~ 400 meV [8 ) and the 
observed FS curvature {t_\_/t\\) is ~ 15% greater in ErTes. 
These differences account for the large intensity observed 
along r — X as one might expect this region to be gapped 
dit Ef based on previous RTes studies. With smaller gaps 
and larger FS curvature, the downward shift of the CDW 
gaps due to the length of the CDW vectors and Px — Pz 
hybridization allows intensity from the top of the CDW 
gaps to be observed at Ef- 

To further illustrate the implications of CDW forma- 
tion on the FS, temperature dependent FS data are 
taken. One significant advantage of ErTes is the fact 
that the room temperature structure is free from any 
CDW Thus this system provides a unique oppor- 

tunity to observe the FS evolution across multiple CDW 
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transitions. The large gaps observed at low temperatures 
for this family allow us to observe the gap opening despite 
thermal broadening effects. Fig. [3^-c show the FS data 
taken at different temperatures for one side of the inner 
FS and Fig. [3]i-g show the temperature dependence of 
Px/Pz bands near the large and small gaps with a CDW 
and without. At T = 200 K only one gap is evident, while 
at T = 300 K the inner FS square is fully closed as ex- 
pected from a structure unmodulated by any CDW [17] . 
To track the temperature evolution of the leading edge 
gap, the spectra were first divided by the temperature 
dependent Fermi-Dirac function convoluted with an en- 
ergy resolution function [21] . Such a procedure allows the 
determination of the center of the energy gap ensuring 
accurate measurements of gap values. In addition, since 
the center of the energy gap is observed below Ep^ CDW 
g'-vectors used in the TB model can be determined di- 
rectly from the data resulting in c* — qcDWi ^ 0.70(0)c* 
and a* — qcDW2 ^ 0.68(5)a*. The CDW wave vectors 
determined from ARPES are in excellent agreement with 
the lattice modulation vectors observed in x-ray data [15] . 
The leading edge gap was then determined by fitting mo- 
mentum distribution curve peaks with Lorentz functions 
and tracking the point of inflection in the fitted band 
dispersion after denoising via wavelet shrinkage. Instead 
of tracking the gaps at /c^^ = and /c^ = 0^ temperature 
dependent data were taken at /c-points where the gap 
maxima with no FS intensity are observed. Fig. [3)i sum- 
marizes the temperature dependent data showing both 
gaps closing. A mean-field order parameter curve scaled 
to the maximum observed gap is also plotted for compar- 
ison. The smaller A2 and larger Ai gaps are observed to 
close at Tc2 ~ 160 K and Tci ^ 280 K, respectively, in 
good agreement with the transport and x-ray data [15] . 
The development of the gaps appears to be second order 
within the experimental uncertainty as no hysteresis has 
been observed. While the closing of the gaps is suggestive 
of a mean- field type behavior, Ai(T) is somewhat sup- 
pressed from the mean-field curve. In addition, it should 
be noted that 2Ai/A:bTci is - 2(2A2/A:bTc2) while the 
area of the FS gapped by Ai is ~ 3 times the area gapped 
by A2. 



The observed g- vectors, the observation of CDW gaps 
below the Fermi-level and the use of an electron-phonon 
coupled model Hamiltonian may suggest the FS plays 
little role in the formation of the CDWs [11 . However, 
caution is advised as the model FS in Fig. [2^ is not the 
one electron eigenvalues resulting from the mixing of the 
different bands. To accurately model the FS, the calcu- 
lated spectral weight {px and Pz eigenvectors) had to be 
used. The FS in Fig. ^ is also a poor match for the 
model FS eigenvalues with no CDWs because the bare 
bands folded back into the reduced 3D Brillouin zone are 
too weak to be observed. Lowering of the ground state 
energy is achieved by gapping the FS and the model spec- 
tral weight suggest the shape of the FS could still play a 
significant role in the CDW formation. 

Both CDWs exist within the same Te plane [15^, thus 
both CDWs modulate the positions of the same Te atoms. 
Hence, ErTes offers a unique opportunity to directly 
study the crossover from quasi- ID to quasi-2D behav- 
ior. Upon initial inspection, each CDW appears uni- 
directional and completely decoupled. However, suppres- 
sion of Ai(T) from the mean- field curve, the discrepancy 
between 2A1//C5TC2 and 2A2//cbTc2 still need to be ex- 
plained. Such discrepancies may arise due to the inter- 
play between the two many body states. Subtle com- 
plexities arising from the crystal structure could inter- 
fere with the delicate balance between the lattice and 
electronic energies, allowing for interactions between the 
two CDWs to arise. More experimental and theoretical 
work is required to explore such possibilities. 
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